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Chapter  1 

Current  ballistic  missile,  weapon  systems  are  capable  of  carrying 
multiple,  independently  targeted  reentry  vehicles  MIRV).  ".’be  effective 
assignment  of  these  independent  weapons  against  individual  targets  is  the 
subject  of  the  research  described  below.  Thxs  research  is  a  continuation  of 
the  work  described  in  SRL-Tli-78-0005 ,  optimal  Targeting  of  Ballistic 
Mi s s i 1 es  in  a  Tiered  Aimpo tnt  System . 

The  earlier  work  described  the  modeling  process  which  transformed  a 
weapons  data  base  and  a  target  data  base  into  a  mathematical  programming 
model  whose  solution  was  an  'optimal'  allocation  of  weapons  to  targets- 
Optimality  was  measured  by  maximizing  expecced  damage  value  initially 
defined  by  a  scoring  process.  Ill  this  paper  we  will  a)  review  the 
previously  reported  work,  b)  describe  the  model  changes  necessary  to 
incorporate  the  concept  of  MIRVed  vehicles,  c)  report  the  computational 
results  of  solving  several  MIRVed  problems,  d)  discuss  potential 
computational  improvements,  and  e)  recommend  avenues  for  implementation. 
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Chapter  2.  Mode]  Formulation  and  Assumptions 

2 . 1  IlnMTRVed  Model  Descriptio n 

There  are  several  assumptions  which  are  made  in  the  transition  from  the 
data  describing  the  basic  target  set  and  the  weapon  arsenal  to  the  mathemat¬ 
ical  programming  model  which  is  used  to  construct  a  weapon  laydown.  To 
appreciate  the  utility  of  the  model  and  understand  its  limitations,  one  must 
relate  these  assumptions  to  the  basic  model  structure.  The  next  sections 
review  these  modeling  assumptions  first  for  the  target  set,  then  for  the 
weapon  array,  and  finally  for  the  structure  of  the  model  which  ties  these 
toge  ther . 

2.2  Target  Set 

The  basic  set  of  targets  are  called  installations  and  are  characterized 
by  the  following: 

a)  location  defined  by  coordinates 

b)  value  or  score  (specified  try  targeting  agency) 

c)  hardness 

d)  damage  requirements 

The  Aimpoint  Generation  Algorithm  (AGA)  described  in  SRL-TR- 7 8-0005 
partitions  these  installations  into  aggregated  targets  called  designeted- 
ground-zeros  (DGZs).  For  a  particular  weapon  type  (tier)  these  DG'is  are 
notiover lapping  target  clusters  meeting  the  criteria  specified  by  the  AGA 
user.  A  DGZ  is  defined  hy  the  following: 

a)  location  (coordinates) 

b)  eligible  weapon  type(s) 

c)  valla-  extracted 


d )  he  i  "h L 


of  burst 
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2.3  Weapon  Data  Set 

A  single  weapon  type  is  assumed  to  be  launched  from  a  uniquely  specified 
launch  site,  at  a  location  defined  by  its  coordinates.  The  assumption  of  a 
single  weapon  type  per  launch  site  is  noL  a  limitation  since  multiple  launch 
sites  can  be  defined  at  tha  same  coordinates.  The  number  of  weapons  at  a 
launch  site  is  determined  by  dei.ir.ing  the  number  of  launch  vehicles  per  site 
and  the  number  of  weapons  per  launch  vehicle.  Any  launch  site  which  has 
more  than  one  weapon  per  launch  vehicle  is  called  a  MIRVed  launch  site. 

Thus,  the  defining  data  elements  for  a  launch  site  are: 

1)  location  (coordinates) 

2)  weapon  type  (W/CEP) 

3)  number  of  launch  vehicles 

4)  number  of  weapons  per  launch  vehicle 

5)  range  limitation  for  different  launch  azimuths  for  the  weapon  type 

2.4  UnMIKVed  Model  _De_fin i trio n 

From  the  target/weapon  data  base  defined  in  the  previous  two  sections,  a 
mathematical  programming  model  can  be  formulated  which,  when  optimized,  will 
solve  the  problem  of  allocating  weapons  to  targets.  The  fundamental  decision 
variable  0^  .)  is  defined  for  all  DC ",  s  which  are  eligible  targets  for 
launch  site  i  under  all  restrictions  such  as  range,  launch  azimuth,  and 
weapon  type.  The  set  of  eligible  DGZs  for  launch  site  i  is  designated  D(i). 
The  decision  variable,  >.  .,  represents  weapon  assignment,  thus: 


1  If  a  Wfinon  from  launch  site  is 
allocated  to  DC.Z  j 

0  it  he  rw  iso 
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The  value  extracted  is  an  output  of  the  AGA  and  represents  the  expected 
collective  damage  to  ail  the  installations  within  a  DGZ  if  one  of  the  desig¬ 
nated  weapon  type(s)  is  detonated  according  to  the  parameters  specified  when 
the  AGA  created  the  DGZ.  The  computation  of  the  value  extracted  assumes  no 
interaction  with  other  DGZs ,  i.e.,  no  collateral  damage.  It  is  an  expected 
value  figure  taking  into  account  weapon  inaccuracies  expressed  as  CEPs  and  a 
user  designated  probability  of  minimum  damage  (P  )  used  in  the  expected 
value  calculation. 

Another  set  of  data  is  required  to  completely  describe  the  characteris¬ 
tics  of  the  target  set  for  model  construction.  Although  the  DGZs  within  a 
tier  or  weapon  class  are  constructed  ar.  nonoverlapping,  the  DGZs  defined  for 
different  tiers  are  either  proper  subsets  of  another  tier  or  are  overlapping 
covers  of  the.  target  set.  This  DGZ/tier  coverage  interaction  is  described 
by  identifying  DGZs  by  tier  and  specifying  for  each  installation  which  DGZs 
include  it.  Since  tier  definition  is  arbitrary,  this  is  a  flexible  method 
adaptable  for  a  wide  range  of  operational  situations. 

To  review  the  definitions  used  to  construct  the  target  set: 

1)  The  DGZ  represents  the  basic  aimpornt  or  target  for  its  associated 
weapon  type. 

2}  The  value-extracted  represents  the  value  of  allocating  the  DGZ 
weapon  type  to  that  DGZ;  no  interaction  with  other  DGZs  or  weapons  occur, 

3)  The  value-extracted  is  an  expected  value  measure  and  incorporates 
probabi  '  iti.es  quant  j  tying  inaccuracies  in  weapon  delivery. 

4)  The  collateral  coverage  between  DGZs  oi  different  tiers  is 
completely  defined  h  ■  including  appropriate  installations  in  the  respective 
DGZs , 

f. 


As  described  previously,  each  DGZ  is  assigned  a  value  which  is  extracted 
by  assignment,  of  a  specified  weapon  system.  Thus,  the  objective  of  maxi¬ 
mizing  extracted  value  (in  a  sense,  maximize  expected  damage)  can  be  mathe¬ 
matically  stated: 

MAXIMIZE  £  S  V.  Aij 

i  i  e  n(i) 

Inherent  in  this  additive  definition  of  the  objective  function  is  the 
assumption  that  at  most  one  weapon  is  assigned  to  each  DGZ.  Multiple  weapon 
assignments  would  require  a  nonlinear  objective  function  since  the  damage  is 
not  additive.  Because  our  environment  is  "target  rich",  i.e.,  there  are 
many  more  targets  than  available  weapons,  the  restiiction  of,  at  most,  one 
weapon  per  DGZ  does  not  markedly  degrade  the  quality  of  the  resulting  weapon 
allocation.  In  addition,  the  definition  of  a  DGZ  is  essentially  a  grouping 
of  installations.  These  installations  are,  in  fact,  the  actual  targets. 

The  DGZ  designations  allow  the  inclusion  of  a  particular  installation  in 
more  than  one  DGZ.  The  nonadditive  characteristic  of  damage  precludes  the 
assignment  of  weapons  to  two  DGZs  which  overlap,  i.e.,  which  contain  some  of 
the  same  installations.  All  of  these  multiple  coverage,  overlap  considera¬ 
tions  can  be  included  in  a  single  type  of  restriction: 

Y  1  *  V  k 

(i,j)  i  Uk) 

where  l(k)  is  the  set  of  launch  site  DGZ  pairs  which  describe  eligible 
weapons  to  a  Dt»Z  containing  installation  k , 

Taken  literally,  this  cover  constraint  would  produce  many  redundant  and 
trivial  constraints.  Thus,  in  implementing  the  model,  only  those  explicit 


non redundant ,  nontrivial  const  taints  are  included. 


The  weapon  supply  is  limited  by  the  number  of  launch  vehicles  (n_.  )  and 
the  number  of  reentry  vehicles  per  launch  vehicle  (lb)  at  each  launch  site. 
For  an  unMIRVed  launch  site,  ib  =  1.  Thus,  the  allocation  for  a  particular 
launch  site  is  limited: 


£  Hi  *  pini  Vi 

i  r  P(1) 

The  unMIRVed  problem  can  be  summarized  as  a  zero/one  integer  programming 
problem : 

MAXIMIZE  £  £  Vi  Ai1 

i  ieD(i) 


subject  to 


(DIP) 


X  Vi  i  i  vk 

(ij)  £I00 


j  e  D  ( i ) 

=  0  or  1 


2.5  Footprinting 

In  the  problem  UIP,  each  reentry  vehicle  is  treated  as  an  independently 
targeted  weapon,  subject  only  to  the  range  and  weapon  type  restrictions 
imposed  in  the  definition  of  the  decision  variables  A{j ■  In  reality,  when 
several  reentry  vehicles  (RV)  are  on  a  single  launch  vehicle,  there  are 
limitations  imposed  on  the.  RVs  as  a  group.  Simply  put,  the  total  energy 
available  for  dispersing  the  group  of  RVs  is  bounded;  implying  that  the 
collective  range  of  the  group  is  limited.  This  limitation  could  be  satisfied 
in  a  variety  oi  ways  such  as  all  but  one  RV  landing  in  a  tightly  clustered 
group  and  the  remain  Lug  RV  using  most  oi  the  energy  lo  hit  a  distant  target 
(Figure  )  .  (n  applying  this  proximity  constraint  to  MTRVed  weapons,  one 
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could  consider  all  possible  sets  of  DGZs  eligible  for  a  specified  MIRV  launch 
vehicle  (MLV),  and  identify  those  which  meet  the  MIRV  proximity  condition. 
This  would  produce  an  unreasonably  large  number  of  target  alternatives  and 
would  be  computationally  difficult  to  implement.  A  way  to  produce  a  rich 
but  manageable  number  of  MIRV  target  alternatives  is  to  pick  a  geometric 
figure  representative  of  a  realistic  MIRV  laydown  which  can  be  varied  by 
simple  parameter  specification,  and  assume  that  any  pattern  of  RV  laydown 
within  fne  figure  is  feasible.  This  figure  is  designated  the  MIRV  footprint . 

An  ellipse  provides  sufficient  flexibility  yet  is  easy  to  implement; 
thus,  it  was  chosen  as  the  MIRV  footprint  mechanism.  A  footprint  ellipse  is 
oriented  witn  its  major  axis  along  the  radius  from  a  LS  and  has  a  dimension 
ratio  which  can  be  set  to  represent  the  desired  pattern  characteristics. 

For  simplicity,  the  footprints  are  all  taken  as  identical,  however,  the 
model  requires  only  that  the  shape  and  dimension  of  footprints  be  identical 
for  all  LVs  deployed  from  the  same  LS .  A  typical  MIRV  footprint  might  have 
an  ellipse  ratio  of  5  with  i  largest  dimension  of  approximately  200  miles. 

Figure  2  shows  the  footprints  for  a  3  weapon  M1KV  ( u  =  3)  targeted 
against  a  small  target  set. 

?.6  Print  Set  and  MIRV  Target  Set 

By  overlaying  a  footprint  on  the  set  of  DGZs  eligib-le  to  a  particular 
1.S,  we  define  those  DGZs  within  the  footprint  as  a  print  set.  Thus,  a  print 
set  is  a  set  of  DGZs  eligible  for  targeting  from  a  single.  MLV  and 
encompassed  within  a  specific  footprint.  It  is  clea'  that  many  print  sets 
with  multiple  redundancies  could  be  defined  by  repositioning  the  footprint 
and  defining  the  set  of  DGZs  covered  as  a  new  print  set.  F.ven  a  scheme  of 
defining  an  arbitrary  grid  and  positioning  the  footprint  on  all  points  of 
tin;  grin  would  generate  an  unmanageably  large  number  ni  print  sets.  The 
technique  used  to  define  print  sets  is  to  position  the  center  ot  a 
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O  SECONDARY  TIER  DGZ'S 


LAUNCH  AZIMUTH 


FIGURE  2.  ILLUSTRATIVE  FOOTPRINTS 


footprint  at  each  eligible  DGZ  and  define  the  covered  DGZs  as  a  print  set. 
Furthermore,  we  will  only  define  print  sets  which  conform  to  the  reciprocity 
condition  given  as: 


(X, 


-  V 


where  x 

and  y  are 

and  S, 

and  R,  are 

k 

k 

downrange  and  crossrange  coordinates 
parameters  of  the  footprint  ellipse, 

Semimajor  axis 


measured  from  LS  k, 

and  R,  and  S,  define: 
k  k 


Semiminor  axis 


Thus,  the  set  of  decision  variables 
set  defined  about  DGZ  i  would  be 


kj 


"i  e  V‘> 


from  launch  site  k  contained  in  the  print 


where  k^(i)  is  the  print  set  defined  about  DGZ  i  for  a  MLV  from  launch 
site  It.  This  method  of  defining  print  sets  provides  a  rich  but 
computationally  manageable  number  of  potential  MIRV  targets.  Figure  3 
illustrates  this  concept  of  a  print  set. 

The  number  of  DGZs  within  a  print  set  depends  on  the  geometry  and 
characteristics  of  a  specific  launch  site/target  set  scenario.  In  general, 
the  number  of  DGZs  in  a  print  set  would  exceed  the  number  of  weapons  on  a 
MLV.  To  distinguish  those  DGZs  specifically  targeted,  the  MIRV  target  set 
will  be  defined  as  that  subset  of  U  DGZs  within  a  print  set  which  the 
weapons  on  a  MLV  are  allocated  against.  This  is  shown  in  Figure  4. 
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FIGURE  3.  FOOTPRINT  PRINT  SET  DGZ  HIERARCHY 
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B.  EXAMPLES  OF  PRINT  SETS 


FIGURE  4.  ILLUSTRATION  OF  A  MIRV  TARGET  SET 
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A.  FOOTPRINT  AND  ASSOCIATED  PRINT  SET 
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2.7  Modeling  the  MIRV  Constraints 

In  characterizing  the.  Ml KV  weapon,  a  hierarchy  of  terms  has  been 
required.  The  footprint  is  a  geometric  figure  whose  size,  shape,  and  orien¬ 
tation  are  dependent  on  the  specific  MLKV  weapon  system  and  launch  site  and 
which  defines  an  area  of  feas'fcle  coverage  by  a  single  MIRV.  A.  print  set  is 
the  specific  set  of  eligible  DGZs  for  a  particular  placement  of  a  footprint. 
A  MIRV  target  set  is  that  subset  of  a  print  set  designated  to  be  hit  by  a 
reentry  vehicle  from  the  MI.V. 

In  formulating  the  MIRV  problem,  the  three  following  conditions  must  be 
satisfied  by  MIRV  target  sets; 

a)  Each  MIRV  target  set  must  be  a  subset  of  at  least  one  print  set. 
Clearly,  a  MIRV  target,  set  can  overlap  two  or  more  print  sets,  and  any  print 
set  can  contain  more  than  one  MIRV  target  set,  but  each  MIRV  target  set  must 
be  a  subset  of  at  least  one  print  set.  (See  Figure  4.) 

b)  Use  of  a  single  weapon  from  a  MI.V  implies  use  of  all  the  weapons 
from  that  MI,V.  That  ls  if  there  are  |i  weapons  on  a  MI.V,  targeting  one 
implies  targeting  IJ  against  a  specific  MI  KV  target  set;. 

c)  MIRV  target  sets  must  be  disjoint;  they  cannot  share  a  curcmot.  I)G 2. 
This  is  implied  by  the  single  coverage  assumption  necessary  to  assure 
linearity  of  the  objective  function. 

These  requirments  will  now  be  translated  into  modifications  oi  the 
•nMTRVed  assignment  problem.  whereas  conditions  (b)  and  (c)  wiii  be  seen  as 
variants  of  supply  and  single  coverage  const  raters,  condition  (a)  has  no 
counterpart  in  the  unMIRVed  formulation.  It  introduces  a  new  constraint. 

For  the  moment,  we  will  drop  Lhe  subset ipl  I  and  consider  a  single  MIKV  l.b 
with  b  weapons  per  I.V. 

I  h 


The  number  of  assignments 


print  set  i.  is 

E 

i  c.?(i) 

Condition  (a)  states  that  an  assignment  is  permitted  to  DGZ  k  only  if  there 

are  at  least  \ i  assignments  in  some  print  set  containing  k,  that  is,  only  if 

y  ’  X .  >  p  for  some  i  f.  P(k) 

I  c.P(i) 

To  distinguish  print  sets  containing  at  least  p  assignments  from  those 
containing  less,  we  define  a  new  variable,  w. ,  by  the  relation 

V  £  Af  “  V  W.  >  0  (la) 

,ir.P(J) 

x^here 


>  0  and  integer 


(lb) 


Then 


'.r  = 

i 


0  if  print  set  i  contains  less  than  h  assignments 
0  if  print  set  i  contains  u  or  more  assignments 


Since  k  (  P(j)  impl  ies  j  f;  P(k),  the  quantity 


E  », 

irP(k) 


is  greater  that  zero  only  if  at  icasL  one  print  set  containing  DGZ  k  contains 
P  or  more  assignments.  The  mathematical  expression  of  (a)  is  therefore 


k 


0 


togeth.  r  with  the  defining  relations  (  I  '  They  relations  (with  the 
subset  i  pt  1  restored  lo  all  quant  il  irs  appearing)  comprise  the  print 


I  h 


constraints.  They  are  consistency  relations  ensui ing  that  MIRV  veapons  are 
allocated  in  proximity  clusters, 

As  indicated,  condition  (c)  arises  from  the  single  coverage  restriction, 
but  adoption  will  require  redefining  the  decision  variables.  Condition  (b) 
will  then  follow  with  a  slight  corresponding  modification  of  the  supply 
constraint  s , 

Condition  (c)  requires  that  MIRV  target  sets  be  disjoint  They  must, 

therefore,  be  identifiable.  To  this  point,  weapons  at  a  single  LS  were 

considered  indistinguishable.  The  decision  variable  A.^  .  was  equal  to  1  if 

any  weapon  from  LS  1  was  allocated  to  ncz  j  and  0  otherwise.  The  single 

coverage  constraint  ensured  that  LS  target  sets  were  disjoint.  If  now  every 

MLV  is  considered  to  be  a  separate  LS  (many  with  the  same  coordinates),  then 

the  corresponding  LS  target  sets  are  MIRV  target  sets  and  will  be  disjoint 

by  the  singlp  coverage  constraint  operating  on  the  expanded  set  o.f:  decision 

variables,  A  where  1  now  denotes  I.S  or  MLV  as  the  case  may  be.  MIRV 
1  j 

target  sets  are  then  identified  by  MLV  and  DGZs  targeted.  And  finally,  this 

modification  enables  the  remaining  requirement,  condition  (b),  to  be 

satisfied  if  the  supply  constraints  ale  written 

n,  i  f  1  I  ,s  an  unMirved  LS 

E>  < 

i  '  I1 1  if  J  is  a  MllU'ed  LS 

for  the  print  constraints  require  at  least  |i  assignments  while  n-a  mor. 
that)  ! 1 1  assignments  are  permitted  by  the  supply  constraints. 

I  / 


The  MIRV  prohlem  can  be  now 

stated : 

MAXIMIZE 

E  E 

i  jED(i) 

■SUBJECT  TO: 

SINGLE  COVER 

2  5  l 

k  -  1, 

(i,j)  e  d(k) 

SUPPLY-UNMTRVED 

E  <  *i- 

I  e  L 

u 

jGD(i) 

SUPPLY-  M1RVED 

E  -  vi 

it-:L 

m 

jcD(i) 

POINT  SET  ELIGIBILITY 

Ewj  -  o 

■icL 

m 

jcR(k) 

k  =  1 j • • • * 

ALLOCATION  WITHIN  A  PRINT  SET 

V  A .  ~V.  W. 

i—*  lk  l  k 

>  0  it:L 

—  m 

kfQ(j) 

k  =  1 

>  0,  integer  A_^  =  0  or  1 


Solution  of  MIP  will  yield  a  laydown  of  weapons  which  triaximizes  the  value 
assigned  to  the  targets  designated  to  be  hit  within  the  supply  constraints 
and  the  MIRV  proximity  constraints  under  the  assumption  that  n'  installation 
will  be  hit  with  more  than  one  weapon. 

The  price  paid  for  incorporating  these  proximity  constraints  is  the 
substantial  increase  in  the  number  of  variables  required  to  represent  the 
MIRV  model.  Table  1  shows  the  magnitude  of  the  increase  for  several  example 
problems.  The  actual  increase  depends  on  launch  site  and  target  geometry 
and  MIRV  characteristics  since  Lhese  determine  the  number  of  print  sets, 


number  of  DGZs  per  print  set,  etc. 


Table  I 


COMPARISON  OF  UNMiRVED  AND  MIRV  PROBLEM  SIZE 


UNMIRVED  MIRVED 


PROBLEM 

SIZE 

SIZE 

1 

97 

X 

530 

1057 

X 

1010 

2 

100 

X 

752 

1 504 

X 

1454 

3 

101 

X 

730 

1561 

X 

1460 

4 

101 

X 

802 

1515 

X 

1465 

5 

99 

X 

684 

136/ 

X 

1318 

6 

99 

X 

666 

13  31 

X 

1282 

7 

9  7 

X 

5  30 

10  57 

X 

1010 
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The  size  of  the  problem  is  difficult  to  estimate  in  general  because  it 
depends  on  the  geometry  of  the  launch  sites  and  targets  which  dictate  target 
eligibility.  The  number  of  constraints  can  be  determined  by: 

])  one  for  each  installation  (single  coverage) 

2)  one  for  each  unMIRVed  launch  site 

3)  one  for  each  MLV 

4)  one  for  each  MIRV  pseudo  launch  site  and  eligible  DGZ  pair 

5)  one  for  each  print  set  and  pseudo  MIRV  launch  site  pair 

Four  (4)  and  five  (5)  imply  a  potentially  large  number  of  constraints  since 
there  is  a  pseudo  launch  site  for  each  MLV  weapon.  For  example,  10  MIVR 
weapons  and  eligible  DGZs  would  yield  at  least  600  constraints. 

Likewise,  the  number  of  variables  increases  substantially  when  the  MIRV 
aspect  is  modeled.  Again  the  actual  number  of  variables  depends  on  the 
geometry  and  weapon  characteristics,  but  is  somewhat  proportional  to  the 
number  of  launch  sites  times  the  number  of  eligible  DGZs.  For  the  unMIRVed 
problem  with  two  launch  sites  and  20  DGZs,  the  number  of  variables  was  about 
40  no  matter  how  many  launch  vehicles  were  at  each  launch  site.  If  earn 
launch  site  has  two  launch  vehicles,  each  with  4  MIRVs,  the  number  of  vari¬ 
ables  is  now  160.  Clearly  the  introduction  of  the  MIRV  restriction,  while 
feasible,  imposes  practical  limitations  on  the  size  of  problems  which  can  be 
solved . 


Chapter  3.  Computational  Results 

3.1  Introduction 

Several  problems  formulated  using  the  models  developed  in  Chapter  2  were 
solved  to  demonstrate  the  utility  ot  the.  models  and  investigate  potential 
practical  computational  bounds.  The  problems  were  derived  from  a  represen¬ 
tative  target/launch  site  data  base  provide.u  by  USAK/SA.  The  data  were 
selected  to  exercise  the  range  of  real  world  attributes  present  in  the  model 

3 . 2  Re  suits 

Tables  2  and  3  provide  detailed  statistics  for  the  solution  ol  two  test 
problems  (TP1,  TP2).  The  two  problems  show  two  important  features:  1)  the 
increase  in  computational  comolexity  when  the  MIRV  aspect  is  modeled,  and  2) 
the  amount  of  value  degradation  when  the  MIRV  is  imposed.  In  TP1  the  size 
of  the  problem  increases  by  a  factor  of  at  least  ten  and  the  solution  time 
by  a  factor  of  over  30  while  the  value  of  the  solution  decreased  by  15%. 

For  Tl'2  the  increase  in  size  was  again  a  factor  of  over  ten  while  the  solu¬ 
tion  time  increased  by  a  factor  of  over  7 .  The  degradation  in  the  solution 
value  was  only  VI. 

Table  4  presents  the  results  for  solving  several  problems  generated  from 
the  representative  data  base.  As  ran  be.  seen,  the  number  of  total  weapons 
has  been  varied,  the  MJRV  configurations  have  been  varied  and  the  launch 
azimuth  has  been  varied.  Experience  with  this  set  of  problems  indicate  that 
the  practical  limitation  on  problem  size  which  can  be  solved  using  Lhe 
Burroughs  in t. eget  / 1  inear  programming  package  TEMPO  is  between  1000  x  1000 

and  1.000  x  .  ■- )  0  0 .  As  expected,  Lhe  . . cut  ia)  i  unease  in  solution  time  vs 

problem  size  dictates  the  maximum  solvable  problem  lor  any  iiiatbema'  leal 
programming  (MP)  package.  The  exar'  limitation  varies  because  of  the  Ml’ 
system  imp! ement at  ion  and  uompul or  maun ( act  nr or. 


Table  2 


TEST  PROBLEM  1 

AREA  -  5°  LAT  BY  8°  LONG 
TARGETS  BY  TIER 

NUMBER  VALUE 

PRIMARY  15  4930 

SECONDARY  21_  4827 

36  9757  TOTALS 


UNMIRVKD  SOLUTION 

WEAPONS  13 

LINEAR  PROGRAM  13  ROWS  BY  26  COLUMNS 

SOLUTION  TIME  .15  MIN  SOLUTION  VALUE  4909 


MIRVED  SOLUTION 

WEAPONS  3x3  4  x  1  =  13 

LINEAR  PROGRAM  177  ROWS  BY  I 61  COLUMNS 


SOLUTION  TIME  (min)  VALUE 

LP  .75  4841.625 

1st  INT.  2.52  3734 

BEST  INT.  4.69  4164 

TOTAL  6* 


♦OPTIMALITY  NOT  VERIFIED 


* 


2  3 


Table  3 


TEST  PRO  HI, EM  2 


AREA 

2°  LAT  BY  20° 

LONG 

TARGETS 

BY  TIER 

NUMBER 

VALUE 

PRIMARY 

50 

19557 

SI 

68 

10461 

S2 

81 

13596 

lv  0 

43614 

TOTALS 

(135  SEPARATE 

TARGETS ) 

UNMT.RVED  SOLUTION 

WEAPONS 

31 

LINEAR 

PROGRAM 

9  7  ROWS 

BY  530  COLUMNS 

SOLUTION  TIME 

.86  MIN 

SOLUTION  VALUE 

MIRVED  SOLUTION 

WEAPONS  7x1  2xb 

4x3  =  31 

LINEAR  PROGRAM  961 

ROWS  BY  915  COLUMNS 

SOLUTION 

TIME  (MIN) 

VALUE 

LP 

6.79 

2010: 

1st;  I NT. 

26.9? 

18702 

BEST  INT. 

37 .  2 

151)6 

TOTAL  140* 


OPT  IMA  I, l  Y  NOT  VEKIEIKO 


A  significant  observation  is  the  "cost  of  MIRVing",  that  is  the  reduction 
it  value  extracted  when  the  MIRV  restrictions  are  imposed.  Over  the  range 
of  problems  solved,  this  cost  essentially  remained  at  about  5%  of  the 
unMIRVed  value  extracted.  This  means  that  if  one  is  able  to  "optimize" 
weapon  allocation,  then  the  cost  of  MIRVing  is  not  too  great. 
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Chapter  4.  Conclusions  and  Recommendations 
The  research  described  in  the  previous  chapters  demonstrates  the  feasi¬ 
bility  of  incorporating  the  limitation  of  MIRV  delivery  systems  into  a 
weapons  allocation  model.  Also,  the  flexibility  and  utility  of  such  a  model 
for  use  in  a  planning  environment  is  indicated  by  the  parameter  variation 
and  associated  computational  results.  The  principal  tradeoff  of  this  major 
model  enhancement  is  the  increased  problem  size  and  resultant  increase  in 
computer  resources  required  to  solve  the  linear/integer  programming  problem. 

No  attempt  was  made  to  research  efficient  solution  algorithms  for  the 
MIRV  problem.  The  standard  linear/integer  programming  package,  TEMPO, 
available  on  the  USAF  Academy  B6700  was  used.  In  fact,  default  options  were 
used  for  setting  algorithm  parameters  so  the  computer  resources  required 
represent  a  worst  case. 

If  the  MIRV  model  were  to  be  used  in  a  planning  or  operational  environ¬ 
ment,  there  are  several  things  which  could  be  done  to  improve  solution  times 
and  to  expand  the  size  of  problems  which  could  be  solved.  First,  there  are 
several  commercially  available  linear/integer  programming  computer  packages 
which  solve  problems  faster  than  TEMPO  (5-10  times)  and  more  efficiently 
store  problems  and  thus  allowing  the  solution  of  larger  problems.  Thus, 
implementing  the  MIRV  model  using  one  of  these  codes  will  increase  the  size 
of  problems  which  can  be  solved  and  improve  solution  times.  Also,  in  similar 
applications,  significant  improvement  is  gained  by  tailoring  an  algorithm  to 
take  advantage  of  the  structure  of  a  specific  problem.  This  could  certainly 
be  the  case  of  the  MIRV  model.  It  is  not  unreasonable  to  presume  that  a 
factor  of  10  or  more  improvement  in  solution  times  and  problem  size  could  be 
gained  with  a  special  purpose  algorithm.  This  would  certainly  be  advisable 
for  routine  use  of  the  model. 
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